Higher Order Local Accuracy by Averaging in the Finite Element Method

Let u^ be a Ritz-Galerkin approximation, corresponding to the solution u of an elliptic boundary value problem, which is based on a uniform subdivision in the interior of the domain. In this paper we show that by "averaging" the values of Uu in the neighborhood of a point x we may (for a wide class of problems) construct an approximation to u(x) which is often a better approximation than UyAx) itself. The "averaging" operator does not depend on the specific elliptic operator involved and is easily constructed.